A deeper look at the colors of the Saturnian irregular 

satellites 



Tommy Grav 
Harvard- Smithsonian Center for Astrophysics, 
MS51, 60 Garden St., Cambridge, MA02138 

and 

Instiute for Astronomy, University of Hawaii, 
2680 Woodlawn Dr., Honolulu, HI86822 
tgrav@cfa. harvard, edu 

James Bauer 

Jet Propulsion Laboratory, California Institute of Technology 
4800 Oak Grove Dr., MS183-501, Pasadena, CA91109 
bauer@scn.jpl. nasa. gov 

February 5, 2008 



Manuscript: 36 pages, with 11 figures and 5 tables. 



1 



Proposed running head: Colors of Saturnian irregular satellites 



Corresponding author: 

Tommy Grav 
MS51, 60 Garden St. 
Cambridge, MA02138 
USA 

Phone: (617) 384-7689 

Fax: (617) 495-7093 

Email: tgrav@cfa.harvard.edu 



2 



Abstract 



We have performed broadband color photometry of the twelve brightest irregular satellites 
of Saturn with the goal of understanding their surface composition, as well as their physical 
relationship. We find that the satellites have a wide variety of different surface colors, from 
the negative spectral slopes of the two retrograde satellites S IX Phoebe (S" = —2.5 ± 0.4) and 
S XXV Mundilfari {S' = -5.0 ± 1.9) to the fairly red slope of S XXII Ijiraq (S' = 19.5 ± 0.9). 
We further find that there exist a correlation between dynamical families and spectral slope, 
with the prograde clusters, the Gallic and Inuit, showing tight clustering in colors among most 
of their members. The retrograde objects are dynamically and physically more dispersed, but 
some internal structure is apparent. 



Keywords: Irregular satellites; Photometry, Satellites, Surfaces; Saturn, Satellites. 
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1 Introduction 



The satellites of Saturn can be divided into two distinct groups, the regular and irregular, based on 
their orbital characteristics. The regular satellites are believed to have been created in a scenario 
similar to the Solar System, accreting out of a circumplanetary disk of gas and dust. They are thus 
on prograde, nearly circular orbits close to the equatorial plane of their host planet. The irregular 
satellites are thought to be captured heliocentric objects, but their origin or capturing processes are 
not well understood. 

The first irregular satellite of Saturn was S IX Phoebe, discovered by W. Pickering in 1898. It 
would not be until a full century later that the next irregular satellite of Saturn was announced, 
when Gladman et al. (2001) reported their discovery of 12 new irregular satellites, ranging in size 



from 5 to 35 km in diameter. Subsequent surveys by Sheppard et al. (2003, 2006) and Jewitt et al. 
(2005) have raised the number of known Saturnian irregular satellites to a total of 35. 



The irregular satellites of Saturn were recognized to be clustering in orbital elements by Gladman 



et al. (2001 ). They suggested that the clustering was due to the breakup of several larger progenitors. 



similar to the creation of the asteroidal families. The Saturnian irregular satellites seem to make 
up two clear prograde families, named after Inuit and Gallic gods respectively, and a much more 
dispersed retrograde group, named after the Norse ice giants. It is not clear whether theses retro- 
grades are fragments of one or more captured progenitors, although the great dispersion of the 18 
known retrograde satellites makes the idea of several progenitors more likely. To test the hypothesis 
that the dynamical families are indeed fragments of larger progenitors Grav et al. { 2003 ) presented 
BVRI photometry of 8 saturnian irregular satellites. The prograde families were indeed found to 
have homogeneous colors, but only 3 of the objects had l-cx uncertainty less than 0.05 across all color, 
making it hard to strongly argue in favor of the hypothesis proposed by Gladman et al. (2001). 

In order to place constraints on the many capturing mechanisms that have been proposed to 
cause the creation of the irregular satellite population, it is necessary to understand the number of 
events needed to create the currently observed dynamical structure. Irregular satellites are generally 
considered to have been captured by their host planets from heliocentric orbits Gladman et al. (2001 ). 
Due to the time reversibility of the pure gravitational three-body capture (sun-planet-satellite) the 
object is only temporarily in orbit around the planet unless some dissipative mechanism exist to 
make the capture permanent. At least five different processes have been proposed: 1) capture 
due to collisions between a temporarily captured object and an existing satellite, or between two 
temporarily captured objects (Colombo and Franklin, 1971; Gladman et al, 2001); 2) capture due to 
a mass increase in the host planet (Heppenheimer and Porco_1977); 3) capture due to gas drag from 



a gas envelope or gas disk around the planet ( Pollack et al. 



1979 



Cuk and Burns 



2004a); 4) capture 



through three-body interaction with other tempor arily captured objects ( Jewitt and Sheppard , 2005); 
5) capture during planetary resonance passage (Cuk and Gladman, 2006). It is evident that a 
better understanding of the physical and dynamical structure of the irregular satellites, as well as 
determining the number of capturing events, is essential in constraining and potentially eliminating 
some of these proposed mechanisms. Note that it is still not clear whether one mechanism dominates 
the capture, or whether the irregular satellites observed today are the result of capture through 
several of the proposed mechanisms. 

The study of the physical characteristics of the irregular satellites is also an integral part in 
determining the origins of the progenitors. Understanding where the irregular satellites originate 
from and how they reached the phase space of the giant planets is essential in order to improve 
our understanding of the transfer of ma tter within the e arly So lar System. Several origins have 



been proposed for the irregular satellites. 



Cuk and Burns 



(2004a) proposed that the Himalia family 
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originated among the Hilda asteroids in the outer main asteroid belt. The recent fly by of S IX Phoebe 
by the Cassini spacecraft just 19 days before orbital insertion into the Saturnian system provided a 
wealth of new information. Johnson and Lunine (2005) argued that the measured density of Phoebe, 
at p ~ 1630 ± 45kgm~^, indicates that Phoebe has its origin in the outer Solar System. Origins 
among the trojan population and local population of planetesimals at the later stages of giant planet 
formation has also been proposed. 

This paper presents optical photometry collected using the GMOS-N imager on the 8m Gemini 
North telescope and the LRIS-R&B imager on the the 10m Keck I telescope in preparation for thermal 
observations of the irregular satellites with the MIPS instrument on the Spitzer Space Telescope. It 
also provides an errata to a measurement of Albiorix published in Grav et al. (2003). 



2 Observations 



Optical photometry of 7 Saturnian irregular satellites in the g', r', i' and z' SDSS-fllters were collected 
through queue observations on November 9th and 10th, 2004 and classical observing on April 13th 
and 15th, 2005, both using the Gemini Multi-Object Spectrograph North (GMOS-N) instrument at 
the 8-m Gemini Observatory on Mauna Kea, Hawaii. GMOS-N has a ~ 5 by ~ 5 arcmin fleld of 
view, with 0.073 arcsec per pixel. All observations with this instrument were done in binned mode 
(0.146 arcsec per pixel), which with seeing in the range 0.75 — 1 arcsec gave excellent sampling of the 
point spread function of our targets. Optical photometry of 8 Saturnian satellites were also collected 
in the Johnson-Kron- Cousins BVRI filters on January 6th, 2005, using the Low Resolution Imaging 
Spectrometer (LRIS) on the 10m Keck I telescope, also located on Mauna Kea, Hawaii. LRIS has two 
detectors, the red and blue side of the instrument, allowing observations in two filters simultaneous. 
The red side has a field of view of ~ 7 by ~ 7 arcmin, with a plate scale of 0.21 arcsec per pixel. 
The blue side has field of view of ~ 9 by ~ 9 arcmin, with a plate scale of 0.135 arcsec per pixels. 



In addition we report the re-reduction of one observation from the Grav et al. (2003) that was found 
to be erroneous. The observational circumstances are given in Table]ll 

The data were collected by observing each object in a RBBRVVRIIR sequence (or the equivalent 
for the SDSS filter set), to make sure that rotational light curve effects did not influence the computed 
colors. All images were taken with sidereal tracking and the exposure time was kept short enough 
to avoid trailing of the satellites. Images in each of the fllters were then combined to increased the 
signal-to-noise ratio of the science targets. 

The collected data was reduced using aperture corrected photometry with the IRAF/DAOPHOT 
package (Tody, 1993 Davis, 1994). Science target aperture radii were selected at ~ 1.3 x fwhm and 
aperture correction values were determined using outer aperture radii of 3 — 4 x fwhm on ~ 8 — 12 
fleld stars. Background flux were calculated using an aperture radius from ~ 4 to ~ 8 times the 
full-width-half- max (fwhm) of the flelds point spread function (psf). Zero point magnitudes, airmass 



and color corrections were determined using a set of Landolt ( 1992 ) standard stars taken on the same 
nights as the science observations. In order to compare the resulting photometry the SDSS flltered 
observations were converted to the BVRI fllter system. Transformation equations between the SDSS 



photometric system {g'r'i'z') and the standard Johnson-Kron- Cousins system (Johnson and Morgan 



1953 Cousins 1978) are given in Fukugita et al. (1996) and Jordi et al. (2006). 
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3 Results 



With this survey we have performed accurate color broadband photometry of the 12 most luminous 
Saturnian irregular satellites. We have observed all but one of the prograde satellites (S/2004 Sll in 
the Gallic family was not observed) and 1/5 of the known retrograde satellites of Saturn. This survey 
thus constitutes the most comprehensive study of color photometry of the irregular satellites to date. 
The results are given in Table |2] and show that the colors and spectral slopes of the Saturnian irregular 
satellites are consistent with C-, D- and P-type asteroids. Solar colors were adopted as B — V = 0.67, 
V-R = 0.36 and V - / = 0.71 ( |Hartmann eTaTj [l982[ |T990| |Hardorp| [l980| [Jewitt and Lalul |l998|. 

The colors were converted to reflectance and the resulting broadband spectra are given in Fig. |4| 
[5]and|6]for the Gallic, Inuit and retrograde clusters, respectively. The spectral slope, 5^ [%perlOOnm], 
of the broadband colors are given by 



2(10 



0.4(C-Cs„„) 



1) 



(1) 



where C and Csun are the colors of the object and sun, respectively. AA is the wavelength difference 
between the two filters that make up the color (Luu and Jewitt, 1996). The value given for 5"i in 



Table |3] are found using the B — I colors of the objects and the sun. Since some of the satellites have 
apparent non-linear features, we also fitted the four reflectances to a linear slope S'2 using a linear 
regression method. Note that for the non-linear slopes of S XXIII Suttungr and S XXX Thrymr 
the fit was done excluding the V filter reflectance. Fitting the S'2 gave similar values to the S'l 
calculation. However, S'2 is generally a better way to derive an average slope solution when the 
normalized reflectance's deviate significantly from a linear solution, as it uses all four points (or in 
some cases three points) rather than the two used in calculating S'l. 

To capture the non-linearity we also fit linear slopes for filter subsets, the BVR filters, S'2{b-r), 
and for the RI filters, S'2{r-i)- All four values for the spectral slope is given in Table |3| Using 



the spectral slope to taxonomy class values given in Dahlgren and Lagerkvist (|1995|), we define 
taxonomy classes for our objects 



(again see Table |3j). It should be noted that the spectral slopes 
used in Dahlgren and Lagerkvist ( |1995 ) are based on the spectral slope in the range 400 — 740nm, 
while our slopes span the range 430 — 900nm. Since most of the slopes are near linear, we can still 
safely use the literature values for comparison. The results are plotted in Fig. [7] and show that the 
Saturnian irregular satellites contain a more or less equal fraction of C-, P- and D-like objects. One 
object, Ijiraq, has a spectrum that is marginally redder than the D-type objects and we define this 
'medium red" -type object (MR-type; 15% < S'2 < 25%). No objects appears to have the 



as an 



ultra-red matter (UR-type, S'2 > 25%) found in the trans-neptunian population (e.q. Jewitt, 2002). 

We also look for evidence for the weak absorption feature near 0.7/im that has been attributed 
to an Fe^"*" Fe^"*". This feature has been shown to appear in spatially resolved spectra (Jarvis 



et al, 2000b) and in ~ 50% of the main-belt C-type asteroids through both reflectance spectra and 
photometry (Vilas 1994; Bus and Binzel, 2002 Carvano et al. 2003). Using data from 27 low-albedo 



asteroids, Vilas (1994) showed that there were a ~ 85% correlation between the 0.7/xm and the 3.0/im 



spectral feature due to water hydration. The weighted mean values of the 0.7/im test for all available 
observations are given in Table |4| 

The observed magnitudes were reduced to absolute magnitudes using the H-G formulation 



H{1, 1, 0) = mv - 5 log(rA) + 2.5 log((l - G)$i(a) + G<l>2(a)) 



(2) 



where the first correction term is the distance modulus, with r and A representing the object's 
distance in AU from the Sun and observer, respectively, and the second correction term is the phase 
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modulus (for a more thorough description see Bowell et al, 1989). Using the available observations 
(both from this paper and literature in general) of the saturnian irregular satellites we fit H{1, 1, 0) 
and G using a Levenberg-Marquardt least square method. These values are given in Table |4] and 
the fits are plotted in Fig. [8] to 10 Note that use of photometry reported to the Minor Planet 
Center (MPC)is generally discouraged for any serious application as it is generally unreliable and no 
information on the actual filters used is available. However, the MFC photometry used in this paper 
to derive phase curve parameters are well known to the authors, as they are all from survey and 



follow-up programs in which the authors are members (Gladman et al. . 2001). These observations 
have a signal to noise ratio of ~ 4 in most cases, thus we assume that the error in the MFC collected 
photometry is 0.25 magnitudes. It should be noted that several of the less observed objects still have 
very poor solutions and additional photometric observations are still badly needed to better constrain 
and utilize the phase curve as a marker for relationships between the irregular satellites of Saturn. 
However, using the fitted value for the absolute magnitude and assuming an albedo of pr = 0.08 
(the average albedo of Fhoebe; Simonelli et al. , 1999 ) we calculate the approximate diameter of 



the saturnian irregular satellites observed in this paper (see Table |4]). It should also be noted that 
without knowledge about the rotational light curve caution should be taken in reading to much into 
the phase curve parameter. For example, data collected by the authors do show that Ymir has a 
large rotational light curve amplitude (~ 0.3 peak-to-peak magnitudes) and, while the phase curve fit 
does not show much scatter around the fit, care is needed when interpreting the result. It should be 
noted that arguments about the inadequacy of the H-G formulation in detecting the opposition surge 
for very small phase angles has been made, and several other groups prefer the use of a combination 



of a linear and exponential function with four free parameters (Belskaya et al., 2006 Kaasalainen 



et al., 2003). With additional observations such models could be used to better model the apparent 



strong opposition surge of Faaliaq and Ymir (see Figure |9] and 10 , but with the sparse photometry 
currently available the 2 free parameter of the H-G formulation remains an adequate method for 
determination of the phase light curve. 

In order to study the coUisional events as a possible mechanism behind the creation of the apparent 
dynamical family structure, we used the Gauss equations to calculate the dispersion velocity needed 



to moved the orbits of the fragments from that of the progenitor (Nesvorny et al., 2003, 2004 Grav 



et al, 2003). We further use the same method as Nesvorny et al. (2004) to compute the size of 



the impactor needed to create such a dispersion velocity assuming that 10% of the impact energy 
is deposited into the fragments as kinetic energy. Using the derived sizes from Table |4] and an 
assumed density of p = 1.6 g/cm^ (Johnson and Lunine, 2005) for all bodies we can determine a 
likely scenario for the collision and fragmentation. Since the collision rates among the satellites or 



with cometary objects are too low to explain the apparent family structures (Nesvorny et al, 2004) 



we use remnant planetesimals from early epochs of the Solar System in an attempt to determine the 



disk mass needed to create the families. We use the same three size distributions as [Nesvorny et al.\ 
(2004), and note that the size distribution based on studies of the in situ formation of the Kuiper 



belt (Stern, 1995 Kenyon and Luu, 1998; Kenyon, 2002) is an unlikely analog to a disk of remnant 



planetesimals at Saturn. Nesvorny et al. (2004) showed that this distribution (called C in their work) 



would catastrophically disrupt even the smallest of the known Saturnian irregular satellites, unless 
the disk have a mass M^^^-^ < lOM®. 



3.1 The Gallic Family 

This dynamical group is centered around 37° inclination and currently four members are known 
(see Table |5]). Their semi- major axes are also tightly clustered, a^^g = 207 — 302 saturn radii (Rs), 
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compared to the large spread among all the known Saturnian satellites, a^vg = 186 — A12Rs. The 
largest member, S XXVI Albiorix, is also the object with the lowest semi-major axis among the 
group, Oj^^g = 271Rs- This provides evidence that any fragmentation must have happened after any 
substantial gas envelope dissipated. All the members have high average eccentricity, e^^^ = 0.46—0.53, 
which together with the tight clustering in inclination, strongly suggest that the four objects have a 
common origin. 

We observed three of the four known members of this family and find, contrary to our earlier 
work (Grav et al. 2003), that the family does not have homogenous colors (see Fig. [TJ. The 



new observations show that Albiorix has a significant steeper spectral slope, at S'2 = 12.5 ± 0.4%, 
compared to S XXI Tarvos and S XXVIII Erriapo, with slopes of S'2 = 5.4±1.0% and S'2 = 5.1±0.7% 
respectively. This results does not support the hypothesis that the three objects are fragments of 
a larger progenitor. However, a closer look at the individual observations show that Albiorix has a 
large span of different spectral slopes, from S'2 = 5.3 ± 0.8% in the observations on January 6th, 
2005, to S'2 = 14.9 ± 0.5% in the observation on April 13th, 2005. It is thus evident that the colors 
of Albiorix varies over its surface. Figure |4] shows the reflectance spectra of these three objects. It is 
easily seen that Tarvos and Erriapo have linear slopes to the la confidence level, while Albiorix has 
a linear slope at a 3 — o" (lower) confidence level in the January 2005 observations. The April 2005 
observation, however, is significantly non-linear with a very steep slope from the B- to the R-filter 
{S'2{B-R) = 19.0 ± 0.7) and a shallower R- to I-filter slope {S'2{r-i) = 5.6 ± 1.2). The two new 
observations of Albiorix reported in this paper are displayed separately in Fig. |4] clearly showing 
their differences. The slope of the January 2005 observation of Albiorix is remarkable similar to the 
slopes of both Tarvos and Erriapo. We thus hypothesize that Tarvos and Erriapo are the largest 
fragments from an impact on Albiorix, leaving a less red crater. The phase curves of this family 
further supports the notion that the objects are indeed of similar origin (see Fig. [s]). Their phase 
curves are all fairly shallow, with a weighted mean phase slope of G = 0.39. The Gallic irregular 
satellites thus span the range for main belt asteroids (G ~ 0.1 — 0.5; Bowell et al. , 1989). Additional 
observations are needed to further constrain the phase curve fits for both Tarvos and Erriapo. 

The question is now whether a cratering event is possible for Gallic satellites. We find that a 
dispersion velocity of ~ 90 m/s needed to displace both the two smaller satellites to their current 
orbits from that of Albiorix. Such an impact could have been the result of ~ 1.25km planetesimal 
hitting with a coUisional velocity of 4.79 km/s. The impact would be just sub-catastrophic, excavating 
a crater with a radius of ~ 12 km. Such a collision has a 50% probability of happening for a disk 
with Mjjg]^ ^ 50 — IOOM0, depending on the size distribution of planetesimals. While this is a 



substantial disk, the many large craters on Phoebe is further evidence of such a massive disk (Porco 



et al. 2005). 



It is important to note that Cuk and Burns (2004b) found that the primordial dynamical dis- 



tribution of this family of objects may have been significantly altered when the Great Inequality 
resonance (^ = 5Xs — 2Aj — 2vu — vjs) swept through the region during the migration of Jupiter and 



Saturn to their prese nt configuration. 



his of course complicates our understanding of this family. 



and as mentioned by |Cuk and Burns ( 2004b[ ) more study of the impact of this resonance on the 



Gallic family is indeed needed. In particular, it would significantly help our understanding of the 
evolution of the irregular satellites to know whether this resonance favors the phase space currently 
inhabited by the Gallic family for capture. If so, it could be that the Gallic cluster is populated 
by objects that during the sweeping of the Great Inequality resonance during planet migration was 
moved into the cluster from a larger region of the surrounding phase space. Such a study, coupled 
with the determination of any variations in color on the surface of Albiorix, would indicate whether 
these objects are indeed fragments from a cratering event on Albiorix, rather than the products of 
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individual captures. 



3.2 The Inuit Family 



The Inuit cluster is centered around ~ 45° inclination and have four known members. Their semi- 
major axes are more spread out, Oavg = 186 — 297Rs , than the Gallic cluster, but is still fairly 
clustered compared to the large spread of all the Saturnian irregular satellites. The eccentricities 
are also clustered, with average values e^vg = 0.29 — 0.35. This points to a common origin for the 
members of this cluster. 

For the Inuit family we observed all the four known members and the broadband colors and 
spectra are shown in Fig. |2]and[5} respectively. The results show that this cluster is very interesting, 
showing both dynamical (through the interaction of its members with different resonances) and 
physical diversity. All the four members of this dynamical family have been reported as possibly 

In order to confirm 



[ Carruba et al. 


2002 


Cuk et al. 


2002 


Nesvorny et al. 


2003 



this we used updated orbits from observations collected by the MFC between 2000-2006 (Marsden, 
personal communications), and numerically integrated the four objects using a Wisdom-Holman 



symplectic mapping scheme (Wisdom and Holman, 1991) for 10^ years. The results confirm that 



S XXII Ijiraq and S XXIV Kiviuq are found to reside in the Kozai resonance, while S XXIX Siarnaq 
resides in the secular resonance g — gs = 0, where g and gs are the satellite and saturnian apsidal 



frequencies. As reported by Nesvorny et al. (2003), we also find that S XX Paaliaq does not seem to 
be associated with the g — gs = secular resonance. 



3.2.1 S XX Paaliaq, S XXIV Kiviuq and S XXIX Siarnaq 

The colors of the three satellites, Siarnaq, Paaliaq and Kiviuq, are very similar with a spectral slope 
range of S'2 ~ 10.0 — 13.2. Paaliaq has a somewhat shallower slope than the two others, but also 
has the largest uncertainties of the three. The colors are very similar, with matching colors values in 
B — V, V — R and V" — / at the 1—, 2— and 3 — a levels, respectively. All three objects also test positive 



for the 0.7/im-feature, based on the method proposed by Vilas (1994). The reflectance spectra of 
these three objects are shown in Fig. [5] and the similarity is easily seen. All the color determinations 
of these three objects given in this paper is found to have linear slopes at a 2a confidence level. 

When comparing the phase curve slope parameter, G, of the Inuit cluster it is seen that both 
Paaliaq and Ijiraq have steeper phase curves than Siarnaq and Kiviuq. The fits for Ijiraq and Kiviuq 
is poor, however, due to only having two accurate photometric measurements each. The phase curve 
parameter derived for Siarnaq, G = 0.27 ± 0.04, is lower than reported in Bauer et al. (2006), but 
the two values are within the la confidence level of each other. The new value puts Siarnaq in the 
mid-range for main belt asteroids {G ~ 0.1 — 0.5; Bowell et a/., 1989), but higher than the usual 
values for C-, P- and D-type asteroids (G ~ 0.09 ±0.09; Harris et al., 1989). The opposition surge of 



0.2 magnitudes seen in Siarnaq by Bauer et al. (2006) is not as narrow in our new fit. In fact none 



of the irregular satellites observed show any strong opposition surge, consistent with a high albedo 



surface (Tholen and Tedesco 1994) 



3.2.2 S XXII Ijiraq 

Ijiraq is distinctly different from the three other members of this dynamical family. Its colors are 
moderately redder than the others and it does not show any indication of the 0.7/im-feature. The 
moderately red color found for Ijiraq does not, to our knowledge, appear anywhere in the inner solar 
system, indicating that its origin is the outer solar system where this color is common. The steep 
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phase curve slope further support the outer solar system as an origin for Ijiraq. With a phase curve 
parameter of G = —0.14 ± 0.22 it is very similar to the phase curves of the Centaur population 
(G 0.13 - 0.04; [Bauer et al] [20031 Izofej [Rousselot efal] [2005|). 



3.3 The Retrograde Satellites 



There are currently 26 known retrograde Saturnian irregular satellites and they show dramatic dy- 
namical and physical variety. The average semi-major axes of these satellites span the range from 
213i?s' to 412Rs and they have average inclinations varying from 138° to 176°. Even the average 
eccentricities show large variety from the low average eccentricity of S XXIII Suttungr, e^^^ = 0.115, 
to the highly eccentric S/2004 S18, e^vg = 0.501. The retrogrades are thus unlikely be the result of 
a single capture and subsequent fragmentation. However, it is possible there exist smaller family 
clusters among the retrograde Saturnian irregular satellites, and using dynamical studies together 
with the photometric colors and phase curves reported in this paper will make searching for such 
families possible. The colors of the 5 retrograde irregular satellites that have been observed in this 
paper are shown in Fig. [3| together with the colors of Phoebe and Ymir (Grav et al. 2003). 



3.3.1 S XIX Ymir, S XXIII Suttungr, S XXX Thrymr 



Comparing the results to earlier observations from Grav et al. (2003) reveals some interesting results. 
Interestingly the B-V color for Ymir has been found to have changed significantly. Furthermore, 
recent observations by the authors to determine the rotational light curve for Ymir show that the 
object has a significant peak-to-peak amplitude of ~ 0.3 magnitudes. It thus seems likely that Ymir 
has significant surface variegation and most likely an irregular shape. The result from [Grav et al. 
(2003) is similar to the colors and slopes of Suttungr and Thrymr found in this paper, both which 



have similar mean orbital elements to Ymir (see table [s]). This leads us to speculate whether Suttungr 
and Thrymr might be fragments broken off Ymir during an impact event. The fact that the two 
fragments have shallow slope (indicating that their surfaces have not been reddened by exposure to 
cosmic rays) might be an indication that the impact may have been fairly recent. 

The broadband colors of Suttungr and Thrymr show a sharp dip in the V filter, as can be seen 
in Fig. |6j The While the slope of Thrymr barely can be fitted to a straight line at the 3a confidence 
level giving S'2 = 1.1 ± 0.8, the dip in the reflectance spectra of Suttungr is significant even at the 
3a confidence level when using the best fit. A similar, but less significant trend is also seen in the 



observation of Ymir reported in Grav et al. (2003). We thus believe that the feature is indeed real 



but confirmation of this feature is needed either through more accurate broadband photometry or 
low-resolution spectra. The LRIS B- and V-filters have a central wavelength of Ac = 436.9 nm and 
Ac = 543.7 nm with fwhm values of AA = 88 nm and AA = 92 nm, respectively. It is therefore 
possible that what is seen is a positive feature at wavelengths just below 400nm that increases the 
flux in the B-filter, rather than an absorption feature in the V filter. Such behavior is known among 
some objects in the cometary population (such as P/Tuttle and IRAS-Araki-Alcock; Tholen and 



Zellner 1984) and the jovian irregular satellites J XI Carme (Luu 1991). Spectra of H2O-, C2O- and 



iVi^a-frost all show a significant upturn in the U-filter (Hapke et al. , 1981 ). All of these types of frost 



have distinctive features in the near-infrared, thus Ymir, Thrymr and Suttungr are prime candidates 
for searching for these features. While near-infrared spectra will be nearly impossible to do on these 
three faint objects, careful preparation might allow for the collection of JHK photometry for these 
two satellites. Photometric observations in the U-filter would also be very helpful to confirm the 
possible upturn at low visible wavelengths. 
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The notion that these smaller bodies are fragments of Ymir is, however, put into question by a 
quick look at their orbital elements. Computing dispersion velocities based on Gauss equations we 
find that Thrymr, S/2006 S4 and S/2004 S8 can be explained with a dispersion velocity of 6V ~ 175 
m/s. Suttungr , however, needs a dispersion velocity of 6V > 300 m/s, due to the large difference in 
eccentricity (from 0.338 to 0.115). These dispersion velocities are again significantly higher than what 



is expected from a fragmentation event (Benz and Asphaug 1999 Pisani et al. 



1999 


Pisani et a/. 


1999 


Michel et al. 



2002). Using the formalism from Nesvorny et al. (2004) and assuming that ~ 10% of the impact 



energy is used to disperse the resulting fragments we find that in order for the four fragments to be 
related to Ymir, the progenitor was impacted by a ~ 1.75km large body. The resulting impact energy 



using the coUisional velocity of 4.63 km/s (Nesvorny et al. 2004) is actually on the same order of 



magnitude as the energy needed to fully disrupt the progenitor. While the dispersion velocity is high 
and needs explaining, the significant amplitude of Ymir and the similarities of the colors supports the 
hypothesis that these objects are indeed the fragments of a catastrophic disruption of a progenitor 
with a radius of ~ 17 km. Again using the Nesvorny et al. (2004) formalism such a collision has a 
probability greater then 50% for a protoplanetary disk masses M^i^k ^ 120 — IGOM^, depending on 
the size distribution of the residual planetesimal disk. 

While we have determined phase curve parameters for Ymir, Suttungr and Thrymr in this paper, 
it should be noted that all fits have problems. Both the fits of Suttungr and Thrymr are very poor, 
with only one accurate photometric observation each (reported in this paper) and a few measurements 
reported to the MPC (for which we have assumed an uncertainty of 0.25 magnitudes). Additional 
photometric observations of these two satellites at a variety of phase angles are necessary in order 
to better constrain their fits, and thereby allowing us to compare the results to the other irregular 
satellites. The phase curve fit for Ymir has a low formal error {G = 0.01 d= 0.06), but caution has 
to be taken. Data collected by the authors, but still in early stages of reduction, show that the 
rotational light curve has a significant peak-to-peak amplitude. It is therefore surprising to find such 
a good fit for the phase curve parameters. However, since the determination of any rotational light 
curve for Ymir is still unconfirmed we stand by the current fit, while cautioning the reader. 

Even with the poor fits it seems clear that Suttungr has a fairly shallow slope [G = 0.68 ± 0.46) 
compared to the steep slope of Thrymr (G = —0.22 ± 0.39). This would suggest that the two 
satellites may have different surface composition and are therefore unlikely to have a common origin. 
Furthermore, the phase curve for Ymir is also steep {G = 0.01 ± 0.06). This together with the high 
dispersion velocity needed to move Suttungr to its present orbit from that of Ymir, we conclude 
that it unlikely that these two objects have a common origin. It should be noted that the difference 
in slope between Thrymr and Suttungr has a significance level of only P = 0.136, making it clear 
that additional observations and study are essential in order to better understand any relationship 
between these objects. 



3.3.2 The neutrally colored objects: S IX Phoebe and S XXV Mundilfari 

Mundilfari is a small moon, d ~ 8km, that has average orbital elements of a = 308-Rs; e^vg = 
0.211 and i^vg = 167.0°. There exists at least two other satellites with similar orbital elements, 
S/2004 S17 (320i?s,0.181, 167.8°) and S/2004 S13 (303i?s,0.259,168.8°) that may have the same 
origin as Mundilfari. It would be extremely interesting to get color photometry of these two faint 
satellites, as a similar slope to that of the larger Mundilfari would be strong evidence of a coUisional 
family within the retrograde irregular satellites. 

Mundilfari is found to have a neutral to blueish spectral slope, S"2 = —5.0 ± 1.9, which is linear 
at the 2a confidence level. It is seen to have a very shallow feature in the V filter, similar to that 
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of Suttungr and Thrymr. Its colors are similar to the colors of Phoebe at the 2a confidence level 
and the slope is linear at the 2a confidence level. The similarity with Phoebe leads us to believe 
that Mundilfari might have Phoebe as its origin. A dispersion velocity of SV ~ 450m/s is needed to 
disperse Mundilfari to its current position, assuming that it is a fragment from a cratering event on 
Phoebe (Grav et a/., 2003 Nesvorny et al.. 2003, 2004). Again assuming that ~ 10% of the impact 



energy is used in the dispersion of fragments and a collision velocity of 5.07 km/s (from Nesvorny 



et al. , 2004 ) the impacting body is 2 km in radius and would excavate a crater with a diameter 



of 24km . Images from the Cassini-Huygens space craft show that Phoebe is covered with craters 



ranging from 100km in diameter and down (Porco et a/., 2005). It should be noted that again the 



dispersion velocity is significantly higher than results from laboratory tests or hydrocode simulations. 

Unfortunately, no other accurate photometric observations exist of Mundilfari. Thus an accurate 
phase curve can not be determined. Using the inaccurate photometric measurements reported with 
astrometry to the Minor Planet Center the phase curve is fitted with an essentially flat slope, but 
the accuracy of the fit makes any comparisons of little worth at this point. Additional observations 
of Mundilfari at both low and high phase angles are necessary to test whether it has the same phase 
curve parameter as Phoebe. 

3.3.3 S XXVII Skathi 

This satellite has a spectral slope of 5"2 = 5.2 ±2.8, which identifies it as a P-type surface. The slope 
is linear and appears featureless at the la uncertainty level. Its colors is very similar to that of Ymir, 
but it's orbital parameters make any association with Ymir unlikely. It should be noted that Skathi 
is currently in the region of phase space where secular resonances might have significantly influenced 
i t, but currently its lon gitude of perihelion, zu, is circulating in the direction of its orbital motion 



Cuk and Burns 



2004b). The relationship to the other irregular satellites in its vicinity, particularly 



S/2004 S19, remains uncertain. Further color photometric observations of the many new retrograde 
satellites are needed to further explore the existence of possible families, which would help constrain 
the number of captures needed to create the current distribution. 



3.4 The lack of ultra red matter among the Saturnian irregular satellites 

It is widely believed that the irregular satellites are captured bodies that were formed in the plan- 
etesimal solar disk independent of the planets. The color distribution of the Saturnian irregular 
satellites give little in suggestions for their source region, but non the less it is generally assumed 
that the irregular satellites (especially those of the outer planets) have an origin in the outer solar 



system, as evident by the density of Phoebe (Johnson and Lunine, 2005). Figure 11 shows the spec- 



tral distributions of possible source populations to the irregular satellites. The data for the main 
belt asteroids and trojans was collected using objects in the SDSS Moving Object Catalog with 



and Trujillo (2006) 



photometry with 5m < 0.05 (Ivezic et al, 2002). Values for the outer Solar System populations were 



computed from observations in Peixinho et al. (2004), Doressoundiram et al. (2005) and Sheppard 



The distributions show that the irregular satellites of Saturn most likely have two sources, one 
in the main asteroid belt, and one in the outer Solar System. The distribution of the Saturnian 
irregular satellites is incompatible with both the Centaur, classic Kuiper belt (CKBOs) and Plutino 



populations, but is in fairly good agreement with the scattered disk population (SDOs; Doressoundi- 



ram et al., 2005). The Centaur population also presents a possible source population if one assumes 
that the it consists of two separate population, one with slopes — 15% per lOOnm and another 
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with slopes 15 — 30% per lOOnm (Peixinho et ai. 2003). A better understanding of the processes 
behind the bimodal color distribution among the Centaurs is needed, however, to confidently present 
its less red population as a source region for the Saturnian irregular satellites. The high number 
of negative sloped, bluish sloped satellites in our sample is hard to explain from the distribution of 
Centaurs, SDOs or any other outer Solar system population, thus leading us to speculate that the 
inner population of Jovian trojans or main belt asteroids supply a fraction of the Saturnian irregular 
satellites. This is of course counteracted by the fact that Phoebe is one of these negative, bluish 
sloped objects. Further study is needed to understand this apparent contradiction. 

However, if we assume that the origin of the Saturnian irregular satellites is in the outer Solar 
Systems it seems unlikely that no object with a surface covered by ultra red matter (URM) would be 
captured. No such ultra red matter is seen among the irregular satellites. Ijiraq, with a spectral slope 
of 5*2 = 19.5 ± 0.9, is the highest sloped Saturnian irregular satellite. This is significant less then the 
URM covered surfaces, S" = 25 — 50% per 100 nm, found among both the Centaur and scattered disk 
populations. The same lack of URM is evident in the cometary nuclei, which are also believed to 
have their origins in the outer solar system (Jewitt, 2002). Using a collisional resurfacing model with 
triggered cometary activity Delsanti et al. (2004) concluded that most Centaurs should have very 
red colors, unless their collisional environment is significantly higher than expected. It is likely that 
the collisional environment of the Saturnian irregular satellites are higher than that of the Centaurs, 
due to the more constrained phase space they cover and the gravitational focusing of Saturn on 
heliocentric objects in its vicinity, thus creating less red surfaces due to increased cratering. A more 
thorough investigation on this issue is beyond the scope of this paper, but is very much needed. 

Since the irregular satellites are beheved to have been captured during the later epochs of giant 
planet formation, it is possible that the lack of URM is simply a result of the time scale of space 
weathering in the outer Solar System being longer than that of giant planet formation. This means 
that as objects were scattered inwards from the trans-Neptunian populations and captured by the 
planets, their surfaces were not yet heavily reddened by space weathering effects. A process is then, 
however, needed to halt further reddening for objects captured by the giant planets. 



Jewitt (2002) hypothesize that the formation of ballistic or rubble mantles by sublimation on 



outer Solar System objects as they evolve into short period comets is a possible explanation for the 
lack of URM for the cometary nuclei. Since the Saturnian irregular satellites never enter the water 



sublimation zone no anticipation of global resurfacing by rubble mantles is expected (Jewitt, 2002) 



Ballistic mantles may form on the surfaces of the Saturnian irregular satellites if enough volatiles, 
in the form of CO2 or other supervolatiles, are exposed near their surfaces. Resurfacing through a 
ballistic mantel of debris, resulting in a more neutrally colored surface, would then be created by 
outgassing from active portions of the satellite surfaces. No active outgassing is, however, known 
among the Saturnian irregular satellites. 

While the low density of Phoebe points to the outer Solar System as the origin, it does not 
necessarily mean that it accreted in the trans-Neptunian populations. It is quite reasonable that the 
low density is due to the object originated in the protoplanetary disk between the planets. While the 
giant planets accreted most of the planetesimals in this region, a significant residual disk could survive 
for some time after the giant planets have built their cores. It is possible that some of these residuals 
could be captured into irregular satellites. No currently known objects is believed to have its origin 
in the region of the giant planets, so it is difficult to determine what the density of such a residual 
planetesimal would be. It seems, however, unlikely that the density would differ greatly from that 
of the more distant trans-Nepunian populations. Sheppard and Trujillo (2006) speculate that the 
irregular satellites, Jovian and Neptunian Trojans and the grey excited hot Kuiper belt objects may 
have a similar origin. Dynamical studies show that these population may be have been populated 
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by objects in the current giant planet domain that were dispersed, transported and trapped in their 
current locations during or just after the migration of the giant planets (Gomes, 2003 Morbidelli 
eralj[2005j). 



3.5 The Norse satellites as the source of lapetus dark side 



lapetus is unique among the satellites of Saturn in that the surface on the hemisphere facing in the 
direction of motion is very dark, on the order of a few percent, while the opposing hemisphere has an 



albedo roughly ten times greater. Cassini data (Buratti et a/., 2005) have shown that while the NIR 



spectra of the brightest terrains are consistent with a nearly three-quarters water-ice composition, 
the darker terrain's 1-5 /xm spectra are consistent with water-ice surface concentrations on the order 
of a few percent at most. Ferric absorption bands have also been seen in the Cassini spectrum of dark 
material. Similar bands were reported earlier in ground-based spectra of the dark side of lapetus, and 
in spectra of the darkest terrain on Hyperion (Vilas et al, 1996 Jarvis et al, 2000a), and were used 



to show a possible common origin for the dark material on the two satellites. The hypothesis of a 
similar origin for the dark material of Hyperion and lapetus was further supported by the successful 
application of a two-component model (composed of icy and dark D-type asteroidal material) by 
Buratti et al. (2002) in the analysis of high-resolution ground-based optical spectra. Hence, it has 
been proposed that the dark material that coats the dark, low-albedo leading side of lapetus might 
have an exogenous origin. Furthermore, its origin may be from among the retrograde irregular 
satellites (here termed the Norse satellites as they are named after the ice giants in norse mythology; 



Soter 1974; Cruikshank et al.. 1983; Buratti and Mosher 



, 1995 


Vilas et al. 


1996; 


Owen et al. 


2001 



and that the same source may possibly coat the surface of Hyperion to a lesser extent. The low- 



albedo surface of lapetus has a spectral slope of Sg ~ 7.2% (Tholen and Zellner, 1983). Phoebe was 



the satellite originally proposed as the common source of the dark material, with Poynting-Robertson 



drag as the potential transport mechanism (Soter, 1974). However, the C-type spectra of this largest 



irregular satellite (Clark et al.. 2005) is not compatible with the theory that Phoebe is the direct 



source of the material (Buratti et al., 2005). Some closer matches exist among the irregular satellites 
from the retrograde group in our color-photometry data. The best candidates for the source of this 
material thus far are Ymir and Skathi, which exhibit spectral slopes of 5-8 percent, but the future 
task remains to compare the material spectroscopically. 



4 Conclusions 

• The Saturnian irregular satellites show a surprising diversity in colors, ranging from neutrally 
colored Phoebe and Mundilfari, to the moderately red colors of Albiorix and Ijiraq. 

• None of the observed satellites exhibit ultrared surfaces, as have been observed in the trans- 
Neptunian populations. 

• The surface of Albiorix appears to be variegated (at a > 3 — cr level). Two of the other Gallic 
family members we sampled have similar spectral slopes as one portion of Albiorix's surface 
terrain, and could possibly be products of a collision with the larger irregular satellite. A 
common origin is furthermore supported by their shared shallow phase curve behavior. 

• We find a strong agreement in the spectral slope of three of the four Inuit family members we 
have sampled, possibly indicative of a common origin. 
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• The surface of Ymir also appears to be variegated and may be the parent body of Thrymr and 
Suttungr, broken apart due to a colhsion with a ~ 2 km body. 

• Our regtrograde sample of satellites show marked diversity. However, two members, Phoebe and 
Mundilfari, both show uniquely blue slopes, while dynamical arguments support the possibility 
that they too may have a common origin. 
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6 Table Captions 



Table [TJ The time, instrument and geometry of the observations reported in this paper are 
shown. The solar, r, and observer. A, distance of the object is given in AU. The phase angle, a, 
is given in degrees. Notes a,b,c are given in order to better understand which observation of each 
object corresponds to values given in the subsequent tables. 

Table [2| The color photometry with la uncertainties are given. Notes a,b,c are given in order to 
better understand which observation of each object corresponds to values given in the other tables. 

Table [3} The spectral slopes of the observations are given. All slopes are given in % per lOOnm 
and all uncertainties are la. Notes a,b,c are given in order to better understand which observation 
of each object corresponds to values given in the other tables. 

Table |4| Physical parameters computed based on the available observations of the irregular 
satellites of Saturn. H{1, 1, 0) and G values of Phoebe taken from Bauer et al. (2006). Mean radius 
of Phoebe is taken from [Porco et al. (2005). 0.7/i-tests are based on formalism from Vilas (1994), 
with Y and N reflecting a positive and negative test, respectively. 

Table |5| : Based on lO^y integrations using a Wisdom-Holman symplectic mapping scheme 
(Wisdom and Holman [1991 ) based the current best fit orbits from the Minor Planet Center (Marsden, 
personal communications) . 
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Table 1: Observational Geometry 



Object 


UT 


r 


A 


a 


Notes 


Albiorix 


LRIS-R&B 2005-Jan-06 10:40 


8.9356 


7.9602 


0.8545 


a 




GMOS-N 2005-Apr-13 05:40 


8.9742 


8.9755 


6.4035 


b 




ALFOSC 2001-Feb-18 21:30 


9.1750 


9.2304 


6.1467 


c, 1) 


Tarvos 


GMOS-N 2004-NOV-09 05:10 


9.0066 


8.6380 


5.9780 




Erriapo 


GMOS-N 2004-NOV-09 04:00 


8.9398 


8.5467 


5.9653 




Siarnaq 


LRIS-R&B 2005-Jan-06 11:47 


8.9435 


7.9686 


0.8796 


a 




GMOS-N 2005-Apr-15 06:30 


8.9867 


9.0185 


6.3844 


b 


Paaliaq 


LRIS-R&B 2005-Jan-06 10:10 


9.0979 


8.1222 


0.8277 




Kiviuq 


GMOS-N 2004-Nov-lO 15:10 


9.0183 


8.6176 


5.8912 


a 




GMOS-N 2005-Apr-15 06:02 


9.0199 


9.0490 


6.3622 


b 


Ijiraq 


GMOS-N 2004-Nov-lO 12:51 


9.0374 


8.6361 


5.8771 




Ymir 


LRIS-R&B 2005-Jan-06 08:10 


9.1545 


8.1811 


0.9353 


a 




GMOS-N 2005-Apr-13 06:04 


9.1871 


9.1817 


6.2573 


b 


Skatlii 


LRIS-R&B 2005-Jan-06 08:50 


9.1157 


8.1428 


0.9602 




Thrymr 


LRIS-R&B 2005-Jan-06 09:10 


8.8941 


7.9190 


0.8776 




Suttungr 


LRIS-R&B 2005-Jan-06 09:36 


8.9335 


7.9585 


0.8789 




Mundilfari 


LRIS-R&B 2005-Jan-06 12:10 


9.0969 


8.1207 


0.8015 
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Table 2: Color Photometry 



Object 


V 


B-V 


V-R 


V-I 


Notes 


Albiorix 


20.64 ±0.01 


0.80 ±0.02 


0.47 ±0.02 


0.76 ±0.03 


a 




21.32 ±0.01 


0.96 ±0.02 


0.52 ±0.01 


0.96 ±0.02 


b 




21.40 ±0.07 


0.79 ±0.12 


0.52 ±0.07 


1.06 ±0.11 


c,l) 


Tarvos 


23.01 ±0.02 


0.78 ± 0.04 


0.43 ± 0.02 


0.82 ±0.03 




Erriapo 


23.32 ±0.01 


0.71 ±0.03 


0.40 ±0.02 


0.86 ±0.02 




Siarnaq 


20.10 ±0.01 


0.87 ±0.01 


0.48 ±0.01 


1.03 ±0.01 


a 




20.72 ±0.01 


0.88 ±0.01 


0.49 ±0.01 


1.02 ±0.01 


b 


Paaliaq 


21.17 ±0.02 


0.86 ±0.03 


0.40 ± 0.03 


0.92 ±0.03 




Kiviuq 


22.83 ±0.01 


0.86 ±0.02 


0.48 ±0.01 


0.98 ±0.02 


a 




22.63 ±0.02 


0.92 ±0.04 


0.50 ±0.03 


0.98 ±0.03 


b 


Ijiraq 


23.40 ±0.01 


1.05 ±0.03 


0.58 ±0.02 


1.09 ±0.02 




Ymir 


21.88 ±0.02 


0.80 ±0.03 


0.45 ±0.02 


0.89 ±0.02 


a 




22.43 ±0.02 


0.77 ±0.04 


0.41 ±0.03 


0.86 ±0.04 


b 


Skathi 


23.80 ±0.05 


0.72 ±0.08 


0.37 ±0.06 


0.88 ±0.08 




Thrymr 


23.70 ±0.04 


0.41 ±0.06 


0.59 ±0.05 


0.86 ±0.07 




Suttungr 


24.03 ±0.04 


0.47 ±0.06 


0.65 ±0.05 


0.78 ±0.06 




Mundilfari 


24.07 ±0.03 


0.58 ±0.05 


0.41 ±0.04 


0.52 ±0.07 
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Table 3: Spectral Slopes 



Object 




S'2 


S'2{B-Ft) 


S'2{R-I) 


Taxonomy 


Notes 


Albiorix 


4.3 ±0.9 


5.3 ±0.8 


11.4 ± 1.5 


-2.2 ± 1.3 


P 


a 




12.7 ±0.6 


14.9 ±0.5 


19.0 ±0.7 


5.6 ± 1.2 


D 


b 




11.1 ±3.7 


12.9 ±4.1 


12.7 ±5.9 


12.3 ±9.5 


D 


c, 1) 






12.5 ±0.4 


17.3 ±0.6 


2.7 ± 1.0 




Mean 


Tarvos 


5.3 ± 1.2 


5.4 ± 1.0 


8.1 ± 1.7 


2.2 ±2.0 


P 




Erriapo 


4.5 ±0.9 


5.1 ±0.7 


4.0 ± 1.7 


6.1 ± 1.6 


P 




Siarnaq 


12.2 ±0.3 


13.2 ±0.3 


14.3 ±0.6 


12.3 ±0.8 


D 


a 




11.8 ±0.8 


12.5 ±0.9 


9.8 ± 1.4 


17.3 ±2.0 


D 


b 






13.0 ±0.3 


13.5 ±0.5 


12.6 ±0.7 




Mean 


Paaliaq 


9.5 ± 1.2 


10.0 ± 1.2 


11.6 ± 1.9 


9.5 ±3.6 


D 




Kiviuq 


10.9 ±0.6 


11.4 ±0.6 


13.3 ±1.0 


8.9 ± 1.2 


D 


a 




12.2 ± 1.1 


12.9 ± 1.2 


17.1 ±2.0 


7.8 ±2.8 


D 


b 






11.8 ±0.6 


14.1 ±0.8 


8.8 ± 1.1 




Mean 


Ijiraq 


17.5 ±0.8 


19.5 ±0.9 


25.4 ± 1.4 


11.2 ±2.0 


MR 




Ymir 


7.4 ±0.8 


7.4 ±0.7 


9.9 ± 1.7 


5.0 ±2.0 


D 


a 




6.0 ±1.3 


6.1 ± 1.2 


6.9 ±2.1 


5.6 ±2.8 


D 


b 






8.1 ±0.5 


10.0 ±0.9 


5.0 ± 1.4 




Mean 


Skathi 


5.3 ±2.7 


5.2 ±2.8 


2.8 ±4.5 


8.9 ±6.0 


P 




Thrymr 


-2.6 ±2.2 


-3.0 ±2.7 






C 




Suttungr 


-3.1 ±2.0 


-3.2 ± 2.2 






C 




Mundilfari 


-6.7 ±2.0 


-5.0 ±1.9 


-1.4 ±3.1 


-11.7 ±4.0 


C 




Phoebe 


-2.6 ±0.3 


-2.5 ±0.4 


-2.5 ±0.7 


-2.8 ±0.8 


c 
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Object 

Albiorix 

Tarvos 

Erriapo 

Siarnaq 

Paaliaq 

Kiviuq 

Ijiraq 



Table 
H(1,1,0) 
10.87 ±0.01 
12.61 ±0.07 
13.27 ±0.15 
10.24 ±0.02 
11.27 ±0.04 - 
12.43 ±0.16 
12.85 ±0.12 • 



4: Physical Parameters 

G Radius (km) 



0.42 ± 0.06 
0.19 ±0.15 
0.57 ±0.34 
0.27 ±0.04 

-0.04 ±0.12 
0.27 ±0.32 

-0.14 ±0.22 



16 
^ 7 
^ 5 
^21 

13 
^ 8 
^ 6 



0.7/im-test 
N 1.060 ±0.005 
N 1.028 ±0.006 

Y 0.974 ±0.011 

Y 0.937 ±0.005 

Y 0.946 ±0.016 

Y 0.973 ± 0.002 
N 1.024 ±0.012 



Ymir 

Skathi 

Thrymr 

Suttungr 

Mundilfari 

Phoebe 



11.81 ±0.02 
14.04 ±0.11 
13.73 ±0.08 
14.08 ±0.08 
14.28 ±0.08 
6.24 ±0.01 



0.01 ±0.06 
0.64 ± 0.52 
-0.22 ±0.30 
0.68 ± 0.46 
0.95 ±0.52 
0.02 ±0.03 



~ 10 
~ 4 
~ 4 
~ 4 
~ 3 
106.6 ± 1.1 



N 0.997 ±0.011 
Y 0.910 ±0.016 
N 1.201 ±0.027 
N 1.337 ±0.041 
N 1.146 ±0.146 
N 1.023 ±0.006 
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Table 5: Mean Orbital Elements 



Object 




^avg 




^avg 


i 


avg 
















^max 










[AU] 








[o] 




[AU] 




[AU] 












[o] 


[o] 


Gallic Satellites 


Albiorix 





1092 





4768 


37 


16 





1075 





1112 





3256 





6432 


28 


26 


45.88 


Tarvos 





1217 





5254 


37 


93 





1189 





1253 





3580 





7067 


27 


48 


47.95 


Erriapo 





1170 





4700 


37 


41 





1148 





1199 





3156 





6405 


28 


46 


45.89 


S/2004 Sll 





1141 





4684 


37 


88 





1118 





1169 





3075 





6469 


28 


80 


46.67 














Inuit SatcUit 


cs 
















Siarnaq 





1195 





2962 


47 


74 





1174 





1220 





0881 





5565 


40 


35 


54.31 


Kiviuq 





0756 





3154 


47 


72 





0753 





0759 





1358 





5543 


39 


27 


53.70 


Ijiraq 





0758 





3047 


48 


08 





0755 





0762 





1091 





5659 


39 


12 


54.03 


Paaliaci 





1002 





3157 


19 


21 





0992 





1015 





1101 





6177 


38 


37 


56.59 












Retrog 


rade Satellites 
















S/2004 S18 





1355 





5006 


138 


42 





1322 





1398 





2576 





7809 


126 


06 


152.71 


Narvi 





1289 





4278 


142 


11 





1266 





1317 





2272 





6673 


132 


40 


152.76 


S/2004 S19 





1226 





3380 


149 


87 





1206 





1251 





1989 





4970 


142 


80 


156.54 


Skathi 





1041 





2737 


151 


93 





1032 





1053 





1792 





3772 


146 


73 


157.03 


S/2006 S2 





1475 





4773 


153 


37 





1432 





1534 





3063 





6599 


143 


48 


161.98 


S/2006 SI 





1251 





1442 


155 


88 





1231 





1273 





0782 





2218 


151 


82 


159.77 


S/2004 S9 





1354 





2442 


156 


21 





1326 





1390 





1408 





3635 


150 


98 


161.01 


S/2006 S3 





1416 





4511 


156 


75 





1378 





1466 





3018 





6080 


148 


39 


163.91 


S/2004 S15 





1289 





1464 


158 


69 





1267 





1315 





0822 





2225 


154 


72 


162.51 


S/2006 S8 





1178 





4642 


159 


23 





1160 





1203 





3526 





5816 


152 


09 


165.70 


S/2004 S12 





1324 





3289 


164 


62 





1298 





1360 





2322 





4325 


159 


70 


169.09 


S/2004 S16 





1494 





1401 


164 


72 





1456 





1540 





0671 





2240 


160 


90 


168.39 


S/2004 S14 





1316 





3709 


165 


04 





1286 





1356 





2625 





4885 


159 


53 


169.78 


S/2004 SIO 





1388 





2546 


166 


53 





1358 





1427 





1679 





3486 


162 


35 


170.34 


S/2006 S6 





1252 





2175 


162 


79 





1232 





1277 





1477 





2960 


158 


76 


166.61 


Mundilfari 





1243 





2112 


167 


01 





1224 





1265 





1527 





2754 


163 


49 


170.48 


S/2006 S5 





1532 





1897 


167 


66 





1490 





1584 





1049 





2863 


163 


79 


171.23 


S/2004 S17 





1293 





1805 


167 


75 





1270 





1322 





1197 





2483 


164 


25 


171.15 


S/2004 S13 





1226 





2592 


168 


82 





1207 





1251 





1901 





3332 


165 


05 


172.43 


S/2004 S8 





1661 





2088 


170 


15 





1602 





1735 





1061 





3223 


166 


31 


173.69 


S/2006 S7 





1519 





4470 


168 


77 





1463 





1590 





3074 





5818 


163 


21 


173.56 


Ymir 





1535 





3376 


172 


86 





1492 





1592 





2311 





4454 


169 


23 


176.48 


S/2006 S4 





1215 





3267 


174 


24 





1196 





1240 





2533 





4013 


171 


09 


177.53 


Phoebe 





0864 





1635 


174 


91 





0860 





0869 





1404 





1883 


172 


52 


177.92 


Suttungr 





1296 





1152 


175 


69 





1276 





1319 





0734 





1612 


173 


34 


178.78 


Thrymr 





1360 





4659 


175 


24 





1328 





1402 





3603 





5672 


171 


97 


178.95 
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7 Figure Captions 



Figure [T| Tlie observ ed colors for the Gallic family. Both the values from this work, as well as 



values from 



Buratti et al. (2005) and Grav et al. (2003) are plotted 



Figure 2} The observed colors for the Inuit family. Both the values from this work, as well as 



values from [Buratti et al (2005) and Grav et al. (2003) are plotted. 

Figure [3j The observed colors for the retrograde Saturnian satellites. Both the values from this 
work, as well as values from Buratti et al. (2005) and Grav et al. (2003) are plotted. 



Figure |4| The computed reflectance from the observations of the Gallic satellites in this paper are 
shown. Icr uncertainties are shown for all observations and the spectra are offset to better compare 
their features. The S'2 (dotted line), S'2(b-r) and S'2{r-i) (both dashed lines) are also shown. 

Figure [5| The computed reflectance from the observations of the Inuit satellites in this paper 
are shown. The fitted values for S'2 are plotted as dotted lines, la uncertainties are shown for all 
observations and the spectra are offset to better compare their features. 

Figure [6| The computed reflectance from the observations of the retrograde satellites. The 



reflectance of Ymir from Grav et al. (2003) is also shown. The dotted lines give the fitted S'2 spectral 



slope, la uncertainties are shown for all observations and the spectra are offset to better compare 
their features. 

Figure [7| The fitted spectral slope, S'2, of the Saturnian irregular satellites observed in this 
paper. Values for Phoebe is calculated from Grav et al. (2003). The slopes of Thrymr and Suttungr 



are fitted by excluding the V-filter observations. 

Figure [8j The available observations are used to determine the lAU phase curve parameters, 
H{1, 1,0) and G, for the three Gallic satellites. The best fit (solid line) and error bars (gray lines) 
are given. 

Figure [9} The available observations are used to determine the lAU phase curve parameters, 
H{1, 1, 0) and G, for the three Inuit satellites. The best fit (solid line) and error bars (gray lines) are 
given. 

Figure |10^ The available observations are used to determine the lAU phase curve parameters, 
H{1, 1, 0) and G, for the five Norse satellites. The best fit (solid line) and error bars (gray lines) are 
given. 

The spectral slope distributions of several possible origin populations for the Satur- 



Figure 11 



nian irregular satellites compared to the distribution found for the irregular satellites in this paper. 
The ultrared material found in the Centaurs and KBOs does not seem to be present in the Saturnian 
irregular satellite population. However, the satellites do generally match the bluer component of the 
Centaurs. The main belt and Jovian Trojan population have been generating using the SDSS Moving 



Object Catalog (Ivezic et al, 2002). The outer Solar System populations have been generated using 



data from Peixinho et al. (2004), Doressoundiram et al. (2005) and Sheppard and Trujillo (2006). 
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Colors of Gallic satellites 



0.6 



■ Albiorix 
• Tarvos 
^ Errlapo 



0.4 



0.2 L 



0.4 



0.6 



0.8 

V-l 



1.0 



1.2 



0.21— 
0.4 



0.6 



0.8 
B-V 



1.0 



1.2 



Figure 1: 



26 



Colors of Inuit satellites 
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Colors of the Norse satellites 
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Reflectance Spectra of the Gallic Satellites 
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Reflectance Spectra of tfie Inuit Satellites 
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Reflectance Spectra of the Norse Satellites 
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Fitted spectral Slopes for the Saturnian Irregular Satellites 
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Phase Curves of the Gallic Satellites 
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Phase Curves of the Inuit Satellites 
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Phase Curves of the Nordic Satellites 
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Spectral slope distribution of different Solar System populations 
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